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ABSTRACT 



Context. Dense molecular cores are studied in order to gain insight into the processes causing clouds to fragment and form stars. In 
this study we concentrate on a region which is assumed to represent an early stage of clustered star formation in a giant molecular 
cloud. 

Aims. The aims of this study are to determine the properties and spatial distribution of dense cores in the relatively quiescent Ori B9 
cloud, and to estimate their ages and dynamical timescales. 

Methods. The cloud was mapped in the 870 yum continuum with APEX/LABOCA, and selected positions were observed in the lines of 
N2H^ and N2D* using IRAM-30m. These were used together with our previous H^D^ observations to derive the degree of deuteration 
and some other chemical characteristics. Moreover, archival far-infrared Spitzer/MIPS maps were combined with the LABOCA map 
to distinguish between prestellar and protostellar cores, and to estimate the evolutionary stages of protostars. 

Results. Twelve dense cores were detected at 870 yum continuum in the Ori B9 cloud. The submm cores constitute ~ 4% of the 
total mass of the Ori B9 region. There is an equal number of prestellar and protostellar cores. Two of the submm sources, which we 
call SMM 3 and SMM 4, are previously unknown Class candidates. There is a high likelyhood that the core masses and mutual 
separations represent the same distributions as observed in other parts of Orion. We found a moderate degree of deuteration in N^H^ 
(0.03 - 0.04). There is, furthermore, evidence for NiH^ depletion in the core SMM 4. These features suggest the cores have reached 
chemical maturity. We derive a relatively high degree of ionization (~ 10"') in the clump associated with IRAS 05405-0117. The 
ambipolar diffusion timescales for two of the cores are ~ 70 - 100 times longer than the free-fall time. 

Conclusions. The distribution and masses of dense cores in Ori B9 are similar to those observed in more active regions in Orion, where 
the statistical core properties have been explained by turbulent fragmentation. The 50/50 proportions of prestellar and protostellar 
cores suggest that duration of the prestellar phase is comparable to the free-fall time. This timescale can be questioned, however, on 
the basis of chemical data on the IRAS 05405-0117 region. A possible explanation is that this survey samples only the densest, i.e., 
dynamically most advanced cores. 

Key words. Stars: formation - ISM: clouds - ISM: molecules - ISM: structure - Radio continuum: ISM - Radio lines: ISM - 
Submillimetre 



1. Introduction 

Most stars form in clusters and smaller groups in the dens- 
est parts of giant molecular clouds (GMCs). The fragmenta- 
tion of molecular clouds results in dense filaments which con- 
tain still denser cores. These cold star-forming cores are best 
detected using far-infrared (FIR) and submillimetre (submm) 
dust continuum. By studying their physical and chemical char- 
acteristics one hopes to understand conditions leading to pro- 
tostellar collapse and the timescale related to this process. 
Furthermore, the distribution and spacing of dense cores place 
constraints on the fragmentation mechanisms (e.g., turbulent 
fragmentation and ambipolar diffusion) and the possible inter- 
action between newly born stars and their surroundings (e.g., 
Megeath et al. 2008) 1. The estimates of the core masses resulting 
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from dust continuum data can also be used to examine the possi- 
ble connection between the mass distribution of dense cores and 
the stellar initial mass function (IMF), a question of great cur- 
rent interest (e.g., Simpson et al. 2008[ ISwift & William s 20081 
IGoodwin et al. 20081 Rathborne et al. 2009. submitted). 

Parameters affecting the cloud dynamics, such as the degree 
of ionization and the abundances of various positive ions are 
chemically related to deuterium fractionation and depletion of 
heavy species (e.g., CO). Besides being important for the core 
dynamics through the magnetic support and molecular line cool- 
ing, these parameters depend on the core history and charac- 
terise its present evolutionary stage. For example, substantial 
CO depletion and deuterium enrichment are supposed to be a 
characteristic of prestellar cores in the pivotal stage before col- 
lapse (Caselli et al. 1999, Bacmann et al. 2002; Lee et al. 2003]. 
Observations (Tafalla et al. 2002. 2004, 2006) and some theo- 
retical models (Bergin & Langer 1997 lAikawa et al. 2005) sug- 
gest, however, that N-containing species such as the chemically 
closely related nitrogen species N2H^ and NH3 (and their deuter- 
ated isotopologues), remain in the gas phase at densities for 
which CO and other C-containing molecules are already de- 
pleted (e.g., IFlower et al. 20051 2006b and references therein). 
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They are therefore considered as useful spectroscopic tracers of 
prestellar cores and the envelopes of protostellar cores. There is, 
however, some evidence that N2H^ finally freezes out at densi- 

cm"^ (e.g., |Bergin et al. 2002 



ties n(U2) - several xlO^ to > 10'' — 



Pagani et al. 2005 2007). Contrary to this, NH3 abundance ap- 
pears to increase toward the centres of e.g., L1498 and L1517B 
jTafalla et al. 2002 2004). Similar result was found by Crapsi et 
al. (2007) using interferometric observations toward LI 544. 

1.1. OhB9 

Most of molecular material in the Orion complex is concentrated 
on the Orion A and B clouds. Star formation in Orion B (L1630) 
takes mainly place in four cl usters, NGC 2023, NGC 2024, NGC 
2068, and NGC 2071 (e.g.. ILada 19921 ILaunhardt et al. 19961 1. 
The Orion B South cloud, which encompasses the star-forming 
regions NGC 2023/2024 is the only site of O and B star forma- 
tion in Orion B (see, e.g.. Nutter & Ward-Thompson 2007| here- 



after NW07). Apart of the above mentioned four regions, only 
single stars or small groups of low- to intermediate mass stars 
are currently forming in Orion B (Launhardt et al. 1996). 

Ori B9 lies in the central part of Orion B and is a rela- 
tively isolated cloud at a projected distance of ~ 40' (5.2 pc 
at 450 pcQ) northeast from the closest star cluster NGC 2024 



(Caselli & Myers 1995) which is the most prominent region of 



current star formation in Orion B. Ori B9 has avoided previ- 
ous (sub)mm mappings which have concentrated on the well- 
known active regions in the northern and southern part of the 
GMC (NW07 and references therein). 

In this paper we present results from the submm continuum 
mapping of the Ori B9 cloud with LABOCA on APEX, and from 
spectral line observations towards three N2H^ peaks found by 
Caselli & Myers (1994) in the clump associated with the low- 
luminosity FIR source IRAS 05405-0117 (see Fig. 2 in CaselH 
& Myers (1994)fl This source has the narrowest CS linewidth 
(0.48 km s"') in the Lada et al. (1991) survey, and the narrowest 
NH3 linewidth (average linewidth is 0.29 km s"') in the survey 
by Harju et al. (1993). A kinetic temperature of 10 K was derived 
from ammonia in this region. 

We have previously detected the HaD^ ion towards two of 
the N2H^ peaks (Harju et al. 20061. These detections suggest 
that 1) the degree of molecular depletion is high and 2) the 
ortho:para ratio of H2 is low, and thus the cores should have 
reached an evolved chemical stage. The high density and low 
temperature may have resulted in CO depletion. This possibility 
is supported by the fact that the clump associated with IRAS 
05405-0117 does not stand out in the CO map of Caselli & 
Myers (1995). 

In the present study we determine the properties and spatial 
distribution of dense cores in the Ori B9 cloud. We also derive 
the degree of deuteration and ionization degree within the clump 
associated with IRAS 05405-01 17. 



The observations and the data reduction procedures are de- 
scribed in Sect. 2. The observational results are presented in 
Sect. 3. In Sect. 4 we describe the methods used to derive the 
physical and chemical properties of the observed sources. In 



' We assume a distance to the Orion star-forming regions of 450 pc. 

^ The N2H+(1 -0) map of Caselli and Myers (1994) shows two sepa- 
rated gas condensations of ~ 0.1 pc in size. The southern condensation 
has a weak subcomponent. The positions of our molecular line obser- 
vations are given in Table 1 of Harju et al. (2006). 



Sect. 5 we discuss the results of our study, and in Sect. 6 we 
summarise our main conclusions. 



2. Observations and data reduction 

2.1. Molecular lines: NiH^ and A/2D+ 

The spectral line observations towards the three above men- 
tioned N2H^ peaks were performed with the IRAM 30 m tele- 
scope on Pico Veleta, Spain, on May 18-20, 2007. The spectra 
were centred at the frequencies of the strongest N2H^(1 - 0) 
and N2D^(2 - 1) hyperfine components. We used the following 
rest frequencies: 93173.777 MHz (NiU+iJFiF = 123 ^ 012), 
Case lli et al._1995 ) and 154217.154 MHz (N2D+(234 ^ 123), 
Ger in et al. 20011 1. Dore et al. (2004), and very recently Pagani 
et al. (2009b), have refined the N2H^ and N2D^ line frequencies. 
The slight differences between the "new" and "old" frequencies 
have, however, no practical effect on to the radial velocities or 
other parameters derived here. The observations were performed 
in the frequency switching mode with the frequency throw set to 
7.9 MHz for the 3 mm lines and 15.8 MHz for the 2 mm lines. 

As the spectral backend we used the VESPA (Versatile 
Spectrometer Assembly) facility autocorrelator which has a 
bandwidth of 20 MHz and a channel width of 10 kHz. The lines 
were observed in two polarisations using the (AB) 100 GHz 
and the (CD) 150 GHz receivers. The horizontal polarisation at 
higher frequency (D150) turned out to be very noisy and was 
thus excluded from the reduction. The channel width used cor- 
responds to 0.032 km s"' and 0.019 km s"' at the observed fre- 
quencies of N2H^(1 -0) and N2D^(2- 1), respectively. The half- 
power beamwidth (HPBW) and the main beam efficiency, tjmb, 
are 26'.'4 and 0.80 at 93 GHz, and 16" and 0.73 at 154 GHz. 

Calibration was achieved by the chopper wheel method. The 
pointing and focus were checked regularly towards Venus and 
several quasars. Pointing accuracy is estimated to be better than 
4 - 6". The single-sideband (SSB) system temperatures were 
~ 150 - 190 K at 93 GHz and ~ 290 - 340 K at 154 GHz. We 
reached an rms sensitivity in antenna temperature units of about 
0.03 K in N2H+(1 - 0) and about 0.05 - 0.07 K in N2D+(2 - 1). 
The observational parameters are listed in Table [1] 

The CLASS programme, which is part of the GAG software 
developed at the IRAM and the Observatoire de GrenobleQ, was 
used for the reductions. Third order polynomial baselines were 
subtracted from the individual N2H^(1 - 0) spectra before and 
after folding them. From each individual N2D^(2 - 1) spectra the 
fourth order polynomial baselines were subtracted before fold- 
ing. Finally, the summed spectra were Hanning smoothed yield- 
ing the velocity resolutions of 0.064 km s"' for N2H+(1 - 0) and 
0.038 km s^' for N2D+(2 - 1). We fitted the fines using the hy- 
perfine structure fitting method of the CLASS programme. This 
method assumes that all the hyperfine components have the same 
excitation temperature and width (Tex and Av, respectively), and 
that their separations and relative line strengths are fixed to the 
values given in Caselli et al. (1995), Gerin et al. (2001), and 
Daniel et al. (2006). Besides Tex and Av this method gives an 
estimate of the total optical depth, Ttot, i.e, the sum of the peak 
optical depths of the hyperfine components. These parameters 
can be used to estimate the column density of the molecule. 



ihttp : //www ■ iram . f r/IRAMFR/GILDAS | 
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Table 1. Observational parameters. 



Molecule Transition 


Frequency 


Instrument 


-fcff 




Resolution 


^throw 




Obs. date 




[GHz] 








["] [kms-'] 


[MHz] 


[K] 




N2H+ y = 1 - 


93.173777^ 


30m/AB100 


0.95 


0.80 


26.4 0.064 


7.9 


150-190 


18-20 May 2007 


N2D+ J = 2-1 


154.21715~ 


30m/CD150 


0.93 


0.73 


16 0.038 


15.8 


290-340 


18-20 May 2007 


Continuum at 870 yum 


APEX/LABOCA 


0.97 


0.73 


18.6 






4 Aug. 2007 



" JlMB = Seff/^cfF, where B^if and F^ti are the beam and forward efficiencies, respectively. 
* Frequency taken from Caselli et al. (1995). 
Frequency taken from Gerin et al. (2001). 



2.2. Submillimetre continuum 

The 870 jum continuum observations toward the Ori B9 cloud 
were carried out on 4 August 2007 with the 295 channel bolome- 
ter aiTay LABOCA (Large APEX Bolometer Camera) on APEX. 
The LABOCA central frequency is about 345 GHz and the band- 
width is about 60 GHz. The HPBW of the telescope is 18'.'6 (0.04 
pc at 450 pc) at the frequency used. The total field of view (FoV) 
for the LABOCA is 11 '.4-. The telescope focus and pointing were 
checked using the planet Mars and the quasar J0423-013. The 
submm zenith opacity was determined using the sky-dip method 
and the values varied from 0.16 to 0.20, with a median value of 
0.18. The uncertainty due to flux calibration is estimated to be 
~ 10%. 

The observations were done using the on-the-fly (OTF) map- 
ping mode with a scanning speed of 3' s"'. A single map con- 
sisted of 200 scans of 30' in length in right ascension and spaced 
by 6" in declination. The area was observed three times, with 
a final sensitivity of about 0.03 Jy beam"' (0.1 M0 beam"' as- 
suming a dust temperature of 10 K). 

The data reduction was performed using the BoA 
(Bolometer Array Analysis Software) software package ac- 
cording to guidelines in the BoA User and Reference Manual 
(2007fl This included flat fielding, flagging bad/dark channels 
and data according to telescope speed and acceleration, correct- 
ing for the atmospheric opacity, division into subscans, baseline 
subtractions and median noise removal, despiking, and filtering 
out the low frequencies of the 1 //-noise. Finally, the three indi- 
vidual maps were coadded. 

2.3. Spitzer/MIPS archival data 

Pipeline (version S 16. 1 .0) reduced "post-BCD (Basic Calibrated 
Data)" Spitzer/MIPS images at 24 and 70 yum were downloaded 
from the Spitzer data archive using the Leopard software pack- 
agfl 

We used the software package MOPEX (MOsaicker and 
Point source EXtractorj^ to perform aperture and point-spread 
function (PSF) fitted photometry on the sources. The point 
sources were extracted using the APEX package (distributed as 
part of MOPEX). 

At 24 fim a 5.31 pixel aperture with a sky annulus between 
8.16 and 13.06 pixels for background subtraction were used. At 
70 jum the pixel aperture was 8.75 pixels and the sky annulus 
ranged from 9.75 to 16.25 pixels. The pixel scale is 2'.'45/pixel 
for 24 jum and 4'.'0/pixel for 70 fj.m. The MIPS resolution is ~ 6" 
and ~ 18" at 24 and 70 yum, respectively. These values corre- 
spond to 0.01 pc and 0.04 pc at the cloud distance of 450 pc. 

fhttp : //www ■ astro . uni -bonn . de/boawiki/Boa | 
^ |http : // ssc ■ spitzer . ca ltech . edu/propkit/ spot/] 
^ |http : //ssc ■ spitzer . c altech . e du/postbcd/mopex .html | 



The aperture correction coefficients used with these settings are 
1.167 and 1.211 at 24 and 70 jim, respectively, as given on the 
Spitzer Science Center (SSC) websitq^l The absolute calibration 
uncertainties are about 4% for 24 jum, and about 10% for 70 //m 
( |Engelbracht et al. 2007[ [Gordon et al. 2007] l. 



3. Observational results 

3.1. Dust emission 

The obtained LABOCA map is presented in Fig. [T] Altogether 
12 compact sources can be identified on this map. A source was 
deemed real if it had a peak flux density > 5cr (i.e., > 0.15 Jy 
beam"') relative to the local background. The coordinates, peak 
and integrated flux densities, deconvolved angular FWHM diam- 
eters, and axis ratios of the detected sources are listed in Table|2] 
The coordinates listed relate to the dust emission peaks. The in- 
tegrated flux densities have been derived by summing pixel by 
pixel the flux density in the source area. The uncertainty on flux 

density is derived from -^cr^^j + cr^, where cr^ai is the uncertainty 

from calibration, i.e., ~ 10% of flux density, and os is the un- 
certainty from flux density determination based on the rms noise 
level near the source area. We have computed the deconvolved 
source angular diameters, ^s, assuming that the brightness distri- 
bution is Gaussian. The values of 0s correspond to the geometric 
mean of the major and minor axes FWHM obtained from two- 
dimensional Gaussian fits to the observed emission which has 
been corrected for the beam size. The uncertainty on has been 
calculated by propagating the uncertainties on the major and mi- 
nor axes FWHM, which are formal errors from the Gaussian fit. 
The axis ratio is defined as the ratio of deconvolved major axis 
FWHM to minor axis FWHM. Both the flux density determi- 
nation and Gaussian fitting to the sources were done using the 
Miriad software package dSault et al. 19951 1. 

Four of the detected sources have IRAS (Infrared 
Astronomical Satellite) point source counterparts, whereas eight 
are new submm sources. We designate the eight new sources as 
SMM 1, SMM 2, etc. The locations of three N2H+(1 - 0) fine 
emission peaks from Caselli & Myers (1994) are indicated on 
the map with plus signs (see Figs.[T]and |3]l. The NaH^ peak Ori 
B9 E which lies 40" east of IRAS 05405-0117 does not coiTe- 
spond to any submm peak (see Fig.|3j. The N2H^ peak Ori B9 N 
lies about 39" southeast of the closest dust continuum peak (see 
Sect. 5.5). One can see that our pointed N2H+/ N2D+ observa- 
tions, made before the LABOCA mapping, missed the strongest 
submm peak SMM 4 located near IRAS 05405-01 17. 



' |http : //ssc ■ spitzer . caltech . edu/mips/ aper corr | 
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Right Ascension (J20QO) 

Fig. 1. LABOCA 870 /jm map of the Ori B9 cloud. The three large plus signs in the centre of the field mark the positions of our 
molecular line observations (see Fig.[3]l. The small plus sign in the upper left shows the dust peak position of IRAS 05412-0105. 
The beam HPBW (18'.'6) is shown in the bottom left. 

Table 2. Submillimetre sources in the Ori B9 cloud. 



Name 


Peak position 
0^2000.0 [h:m:s] (52ooo.o [° 


:':"] 


^ peak 
870 

[ly beam" ] 


S 870 

[ly] 


8, 
["] 




Axis ratio 


IRAS 05399-0121 


05 


42 


27.4 


-01 


19 


50 





81 


2.7 


±0.3 


30 ± 


4 


1 


3 


SMM 1 


05 


42 


30.5 


-01 


20 


45 





41 


3.0 


±0.3 


57 ± 


7 


2 


5 


SMM2 


05 


42 


32.9 


-01 


25 


28 





21 


0.7 


±0.1 


26 ± 


5 


3 


8 


SMM 3 


05 


42 


44.4 


-01 


16 


03 


1 


14 


2.5 


±0.4 


19 ± 


3 


1 


5 


IRAS 05405-0117 


05 


43 


02.7 


-01 


16 


21 





19 


0.9: 


to.r 


45 ± 


8" 


1 


2 


SMM 4 


05 


43 


03.9 


-01 


15 


44 





28 


1.3 


±o.r 


34 ± 


6 


2 


3 


SMM 5 


05 


43 


04.5 


-01 


17 


06 





16 


0.6 


±0.1 


38 ± 


4 


1 


8 


SMM 6 


05 


43 


05.1 


-01 


18 


38 





26 


2.5 


±0.3 


92 ± 


35 


4 


1 


Ori B9 N 


05 


43 


05.7 


-01 


14 


41 





16 


1.0 


±0.1 


47 ± 


5 


1 


2 


SMM 7 


05 


43 


22.1 


-01 


13 


46 





32 


0.8 


±0.1 


24 ± 


4 


1 


8 


IRAS 05412-0105 


05 


43 


46.4 


-01 


04 


30 





17 


0.5 


±0.1 










IRAS 05413-0104 


05 


43 


51.3 


-01 


02 


50 





66 


0.9 


±0.2 


25 ± 


8 


1 


6 



" These values include both the IRAS 05405-01 17 and Ori B9 E (see text and Fig.©. 
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3.2. Spitzer/MIPS images 



The retrieved Spitzer/MIPS images at 24 and 70 fim are pre- 
sented in Fig. |2] All four IRAS sources that were detected by 
LABOCA are visible at both 24 and 70 fim (IRAS 05412-0105 
and 05413-0104 northeast from the central region are outside the 
regions shown in Fig. From the new submm sources, SMM 
3 and SMM 4 are also visible at both 24 and 70 jum, while there 
is 24 fim source near SMM 5 which is not detected at 70 fim. 
The rest of the submm sources are visible at neither of the wave- 
lenghts. 

In Table [3] we list the sources detected at both 24 and 70 
fim. In this table we give the 24 and 70 yum peak positions of 
the sources and their flux densities at both wavelengths obtained 
from the aperture photometry. The Icr uncertainties on flux den- 
sities were derived as described in Sect. 3.1, i.e. as a quadratic 
sum of the calibration and photometric uncertainties. 

3.3. NxH^ andN2[y 

The three positions of our molecular line observations are indi- 
cated in Figs. [1] and [3] These positions correspond to the three 
N2H+(1 - 0) peaks found by Casefli & Myers (1994) (see also 
Sect. 1.1). The Hanning smoothed N2H+(1 - 0) and N2D+(2 - 1) 
spectra are shown in Figs. |4] and |5] respectively. The seven hy- 
perfine components of N2H^(1 - 0) are clearly resolved towards 
all three positions. The N2H^(1 - 0) spectra towards Ori B9 E 
and Ori B9 N show additional lines, which can be explained by 
N2H^(1 - 0) emission originating at a different radial velocity 
(see Fig. gll. In the case of Ori B9 E, the additional N2H+(1 - 0) 
lines originate in gas at a radial velocity of 1.3 km s"', whereas 
towards Ori B9 N the additional gas component has a Vlsr of 
2.2 km s"'. These velocities are ~ 7 - 8 km s"' lower than the 
average velocity of the Ori B9 cloud, suggesting that they are 
produced by a totally different gas component. We checked that 
the additional components are not caused by, e.g., a phase-lock 
failure by summing randomly selected subsets of the spectra. All 
sums constructed in this manner showed the same features with 
equal intensity ratios. The "absorption"-like feature at ~ 20 km 
s ' in the N2H^(1 -0) spectrum of Ori B9 N is an arfefact caused 
by the frequency switching folding process. 

Only the strongest hyperfine group of N2D^(2 - 1) was de- 
tected. The relatively poor S/N ratio hampers the hyperfine com- 
ponent fitting. Towards Ori B9 N, the additional velocity com- 
ponent at ~ 2.2 km s ' was also detected in N2D"^(2 - 1). 

In Table a we give the N2H+(1 - 0) and N2D+(2 - 1) fine 
parameters derived from Hanning smoothed spectra. The LSR 
velocities and linewidths (FWHM) are listed in columns (2) and 
(3), respectively. The total optical depth and excitation tempera- 
ture for the lines are given in columns (6) and (7), respectively. 
The excitation temperatures, Tex, of the N2H^(1 - 0) transition 
were derived from the antenna equation 



IRAS 05405-0117 



If 



43'"00^ 42"45= 42"'30= 42" 15^ 

Right Ascension (J2000) 




43°00= 42"45= 42"30' 

Right Ascension (J2000) 



4gmjgs 



Fig. 2. Spitzer/MIPS 24 jum (top) and 70 yum (bottom) images of 
the central part of the Ori B9 cloud. The logarithmic colour scale 
range from 33.4 to 729.5 MJy sr"' and 27.9 to 1357.4 MJy sr^' 
in the 24 and 70 fim images, respectively. 



(1) 



where 77 is the beam-source coupling efficiency, h is the Planck 
constant, v is the transition frequency, is the Boltzmann con- 
stant, Tbg - 2.725 K is the cosmic microwave background 
(CMB) temperature, and the function F(T) is defined by F(T) = 
^ghv/keT _ j -y^g assumed that 77 = jjmb, and we used the main 

beam brightness temperature, Tmb - ^mb'^A' '■^^ optical 
thickness, r, of the brightest hyperfine component. The uncer- 



tainty of Tex has been calculated by propagating the uncertainties 
on Tmb and t. 

The total optical depth of the N2D^(2 - 1) line cannot be 
calculated directly because the hyperfine components are not re- 
solved in the spectra. We estimated the total optical depth in the 
following manner First, we calculated the optical depth of the 
main hyperfine group of N2D^(2 - 1) from the antenna equation 
using the Tmb obtained from a Gaussian fit to the group of 4 
strongest hyperfines. In the calculation, we adopted the excita- 
tion temperature of the N2H^(1 - 0) lines. Second, the total op- 
tical depths of N2D^(2 - 1) were calculated taking into account 
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Table 3. Spitzer 24/70 //m sources in Ori B9. 





24 jjm peak position 


70 yum peak position 


S 24 


S 70 


Name 


Q'2000.0 [h:m:s] 


^2000.0 1° '■''■"] 


Q'2000.0 [h:m:s] 


(52000.0 [°:':"] 


[Jy] 


[Jy] 


IRAS 05399-0121 


05 42 27.6 


-01 20 01 


05 42 27.7 


-01 19 57 


1.3 ±0.05 


24.4 ± 2.4 


SMM3 


05 42 45.3 


-01 16 14 


05 42 45.1 


-01 16 13 


0.005 ± 0.0002 


3.6 ±04 


IRAS 05405-0117 


05 43 03.1 


-01 16 29 


05 43 03.0 


-01 16 30 


1.3 ±0.05 


2.7 ±0.3 


SMM4 


05 43 05.7 


-01 15 55 


05 43 05.6 


-01 15 52 


0.036 ± 0.001 


1.3 ±0.1 


IRAS 05412-0105 


05 43 46.3 


-01 04 44 


05 43 46.1 


-01 04 43 


0.6 ± 0.02 


2.5 ±0.3 


IRAS 05413-0104 


05 43 514 


-01 02 53 


05 43 51.3 


-01 02 51 


0.2 ±0.01 


17.9 ± 1.8 



that the main group correspond to 54.3% of Ttot- The uncertainty 
on Ttot has been calculated by propagating the uncertainties on 
Tub and T^^. 



Table 5. IRAS flux densities in Jy. 




Right Ascension (J2000) 

Fig. 3. Blow-up of Fig. [U showing the IRAS 05405-0117 clump 
region. The large plus signs indicate the positions of our H2D^, 
NaH^, and N2D+ observations towards three condensations 
shown in Fig. 2 of Caselli & Myers (1994). Also shown are the 
24 fim peak positions of SMM 4 and IRAS 05405-01 17, and the 
24 yum peak near SMM 5 (small green plus signs, cf. Fig.|2]). The 
beam HPBW (18'.'6) is shown in the bottom left. 



4. Physical and chemical parameters of the sources 

4.1. Spectral energy distributions 

The 24 and 70 fim flux densities together with the integrated flux 
densities at 870 fim were used to fit the spectral energy distribu- 
tion (SED) of SMM 3 and SMM 4. For the IRAS sources also 
the archival IRAS data were included. The flux densities in the 
12, 25, 60, and 100 fim IRAS bands are listed in Table |5] The 
derived SEDs for SMM 3, SMM 4, and IRAS 05405-0117 are 
shown in Fig.|6l 

For all six sources detected at three or more wavelengths, 
the data were fitted with a two-temperature composite model. 
The parameters resulting from the fitting are given in Table |6] 
We adopted a gas-to-dust mass ratio of 100 and dust opacities 
corresponding to a MRN size distribution with thick ice mantles 
at a gas density of «h = 10^ cm"^ ( [Ossenkopf & Henning 1994) l. 



Name 


5,2 


525 


S 60 


S 100 


IRAS 05399-0121 


0.25 


1.59 


22.94 


45.93 


IRAS 05405-0117 


0.40 


1.55 


3.75 


19.67 


IRAS 05412-0105 


0.26 


0.65 


1.66 


73.58 


IRAS 05413-0104 


0.25 


0.31 


17.33 


59.46 



The total (cold-t-warm) mass and the integrated bolometric lumi- 
nosity are given in columns (2) and (3) of Table|6] respectively. 
The temperatures of the two components are listed in columns 
(4) and (5). In columns (6) and (7) we give the mass and lumi- 
nosity fractions of the cold component vs. the total mass and lu- 
minosity, and in column (8) we list the ratio of submm luminos- 
ity (numerically integrated longward of 350 //m) to total bolo- 
metric luminosity (Lsubmm/i'boi)- Column (9) list the normalised 
envelope mass, Mtot/iJ^^p which is an evolutionary indicator in 
the sense that it correlates with the protostellar outflow strength 
(i.e., with the mass accretion rate), and thus decreases with time 



(Bontemps et al. 1996 1. In column (10) we give the source SED 
classification (see Sect. 5.1). In all cases the mass of the warm 
component is negligible (~ 10"^ - lO"'* M0) and thus the bulk 
of the material is cold (Mcoid/A^tot ~ !)■ 

4.2. Linear sizes, mass estimates and densities 

The linear sizes (radii R = 9sd/2) were computed from the an- 
gular FWHM sizes listed in Table |2] 

The masses of the cores (gas+dust mass, Mcont) are calcu- 
lated from their integrated 870 /im continuum flux density, S^ju, 
assuming that the thermal dust emission is optically thin: 



S 870^^ 



(2) 



where d is the distance, and B87o(7^d) is the Planck function with 
dust temperature T^. For all IRAS sources as well as for SMM 3 
and SMM 4 we adopted the dust temperatures resulting from the 
SED fitting (see Table |6l column (4)). For all the other sources 
it was assumed that - 10 K. The assumed dust tempera- 
ture of 10 K is justified by the estimates obtained from NH3 
( [Harju et al. 1993 1 and is commonly adopted for starless cores. 
The assumption that = Tkin, where Jkin is the gas kinetic 
temperature, is likely to be valid at densities n(H2) > 10^ cm"^ 
(Burke & Hollenba ch 19831 1. The opacity per unit mass column 
density at A = 870 fim is assumed to be /cgvo ^ 0.17 m^ kg"'. 
This value is interpolated from Ossenkopf & Henning (1994, see 
Sect. 4.1). The value 1/100 is adopted for the dust-to-gas mass 
ratio, R^. 

The virial masses, Mvh, of IRAS 05405-0117 and Ori B9 
N have been estimated by approximating the mass distribution 
by a homogenous, isothermal sphere without magnetic support 
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Table 4. N2H^(1 - 0) and N2D^(2 - 1) line parameters derived from Hanning smoothed spectra. The integrated hne intensity 
(J r^(v)dv) includes all the hyperfine components in the case of N2H^, whereas for N2D^ only the main group is included. 



Line/Position 


vlsr [km s '] 


Av [km s ' ] 


n 


[K] 


J rX(v)dv [Kkms-'] 




T 

^ ex 


[K] 


N2H+(1 - 0) 


















IRAS 05405-0117^ 


9.228 ± 0.001 


0.290 ± 0.002 


2.37 


±0.04 


3.97 ± 0.05 


6.1 ±0.03 


6.8 ± 


0.07 


Ori B9 E* 


9.163 ±0.002 


0.298 ± 0.005 


1.46 


+ 0.04 


2.26 ± 0.03 


3.5 ± 0.5 


6.1 : 


tO.3 


Ori B9 n3 


9.149 ±0.003 


0.261 ± 0.008 


1.82 


±0.09 


2.29 ± 0.04 


2.2 ±0.8 


8.3: 


t 1.4 


N2D+(2-~) 


















IRAS 05405-0117 


9.414 ±0.012 


0.319 ±0.027 


0.31 


±0.04 


0.13 ±0.01 


0.26 ±0.01^ 


6.8 ± 


0.07 


Ori B9 E 


9.285 ±0.013 


0.194 ±0.029 


0.27 


±0.03 


0.09 ± 0.01 


0.28 ± 0.03~ 


6.1 : 


tO.3 


Ori B9 N" 


9.255 ± 0.009 


0.136 ±0.021 


0.28 


± 0.07 


0.09 ±0.01 


0.16 ±0.05'' 


8.3: 


t 1.4 



" Caselli & Myers (1994) derived Vlsr = 9.209 ± 0.003 km s"', Av = 0.313 ± 0.008 km s"', and t,„, = 4.594 ± 0.825. 
* For the other velocity component hyperfine fit yields vlsr = 1.310 ± 0.013 km s"', Av = 0.436 ± 0.035 km s"', and Ttot = 6.6 ± 2.0. 
For the other velocity component vlsr = 2.219 ± 0.006 km s"', Av = 0.378 ± 0.021 km s"', T'^ = 0.52 ± 0.03 K, Tm = 1.0, and 
Tex = 5.9 ± 0.17 K. 

Ttot calculated by taking into account that the main hyperfine group corresponds to 54.3% of the total optical depth. Tn,ain group is calculated 
using Tmb from Gaussian fit to the main group and from N2H+(1 - 0). 

' For the other velocity component vlsr = 2.290 ± 0.044 km s"', Av = 0.187 ± 0.090 km s"', T* = 0.22 ± 0.06 K, and T,ot = 0.12 ± 0.01 
(Ttot is calculated as described in footnote d). 



Table 6. Results of the SED fits. 





M,ot 




^cold 


T 

warm 






Mto,/Lp 




Source 


[Mq] 


[Lq] 


[K] 


[K] M,m/M,„, 


-^^cold/^bol 


^submm / ^bol 


[M0/L^6] 


Class 


IRAS 05399-0121 


2.8 ±0.3 


21 ± 1.2 


18.5 ±0.1 


103.9 ±0.2 ~1 


0.90 


0.02 


0.45 


o/i 


SMM3 


7.2 ±2.2 


3.5 ±0.2 


11.8±0.9 


36.9 ±0.2 ~ 1 


0.74 


0.11 


5.57 





IRAS 05405-0117 


1.6 ±0.2 


6.4 ± 0.4 


16.1 ±0.1 


112.4 ±0.4 ~1 


0.69 


0.03 


0.53 





SMM4 


3.8 ±0.2 


1.7 ±0.2 


11.6 ±0.2 


51.4 ±4.9 ~1 


0.76 


0.11 


2.76 





IRAS 05412-0105 


1.3 ±0.3 


5.8 ±0.6 


17.0 ±0.3 


127.8 ±0.5 ~1 


0.86 


0.03 


0.45 





IRAS 05413-0104 


1.0 ±0.2 


13.2 ± 1.2 


20.3 ± 0.4 


152.3 ±1.1 ~1 


0.97 


0.02 


0.21 






and external pressure (see, e.g., Eqs. (1) and (2) in Chen et 
al. (2008) where the linewidth of N2H^ is used). The result- 
ing virial masses are about 4.3 Mq for IRAS 05405-0117 and 
2.8 Mq for Ori B9 N. The corresponding Mcont/A^vir ratios are 
about 0.3 and 1.4. Note that it is usual that protostellar cores, 
like IRAS 05405-0117, appears to be below the self-gravitating 
limit (Mcont/A^vir = 0.5), though they are forming stars (e.g., 
lEnoch et al. 20081 1. Since Ori B9 N appears to be gravitation- 
ally bound, it is probably prestellar. There are several factors 
that would lead to virial masses being overestimated. For ex- 
ample, using the radial density profile with power-law indices 
p - I - 1.5 (see Sect. 5.3), would lead to Myh being reduced by 
factors of 1.1 - 1.25. 

The volume-averaged H2 number densities, (n(H2)), were 
calculated assuming a spherical geometry for the sources and 
using masses. Meant, and radii, R, estimated for them from the 
dust continuum map. The obtained radii, masses, and volume- 
averaged H2 number densities are given in colums (2), (3), and 
(5) of Table |7] respectively. 

4.2.1 . Total mass of tPie Ori B9 region 

We made an estimate of the total mass in the region by us- 
ing the near-infrared extinction mapping technique (NICER, 
ILombardi & Alves 20011 ). In this technique, the near-infrared 
colors, namely H - K and J - H, of the stars shining through 
the dust cloud are compared to the colours of stars in a nearby 
field that is free from dust. The reddened colours of the stars 
behind the dust cloud can then be interpreted in terms of extinc- 
tion due to the relatively well-known ratios of optical depths at 
JHK wavelengths (for further details of the method, we refer to 



Lombardi & Alves (2001)). To implement the method, we re- 
trieved JHK photometric data from the 2MASS archive, cover- 
ing a 30'xl9' region centred at {a, 6)s2m) = (5:43:00, -01:16:20). 
Applying NICER to these data yielded an extinction map with 
the resolution of FWHM=2'5, indicating extinction values of 
Ay = 8 ... 12 mag at the positions of the detected sources. 

The total mass of the region was calculated from the de- 
rived extinction map by summing up the extinction values of 
all pixels assuming the gas-to-dust ratio of A^(H) =2-10^' cm"^ 
mag"' (IBohlinetal. 1978l l. and the mean molecular weight per 
H2 molecule of 2.8. The total mass of the region resulting from 
the calculation is 1400 M0. The total mass of the cores within 
the region implied by the submm dust emission data is only ~ 50 
M0, about 3.6% of the total mass in the region. 

4.3. Column densities, fractional abundances, and the 
degree of deuterium fractionation 

The H2 column densities, A^(H2), towards the submm peaks and 
the positions selected for the line observations were calculated 
using the following equation: 

y^dust 

A^(H2) = — . (3) 

/g^Q is the observed dust peak surface brightness, which is re- 
lated to the peak flux density via 1 Jy/18'.'6 beam = 1 .085 ■ 10"'^ 
W m"^ Hz"' sr"'. yUHj = 2.8 is the mean molecular weight per 
H2 molecule, and nin is the mass of the hydrogen atom. The 
same dust temperature values were used as in the mass estimates 

(Eq. my 
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Table 7. Linear radii, masses, and H2 column and volume-averaged number densities of all detected submm sources. 





R 




A'(H2) 


<"(H2)> 


Source 


[pc] 


[Mq] 


[10^- cm-2] 


[10' cm-3] 


IRAS 05399-0121 


0.03 


3.6 


2.61 


7.5 


SMM 1 


0.06 


11.8 


3.89 


3.1 


SMM2 


0.03 


2.8 


1.99 


5.8 


SMM 3 


0.02 


7.2 


7.80 


51.1 


IRAS 05405-0117 


0.05^ 


1.5" 


0.76 


0.7° 


SMM 4 


o.o" 


3.8 


1.98 


3.4 


SMM 5 


0.04 


2.3 


1.51 


2.1 


SMM 6 


0.10 


9.9 


2.45 


0.6 


Ori B9 N 


0.05 


3.9 


1.51 


1.7 


SMM 7 


0.03 


3.1 


3.01 


6.5 


IRAS 05412-0105 




0.8 


0.63 


_b 


IRAS 05413-0104 


0.03 


1.0 


1.85 


2.2 



" These values include both the IRAS 05405-01 17 and Ori B9 E. 

* Deconvolving the angular size was not possible, and thus the radius and number density could not be estimated. 



The N2H^ column densities were calculated using the equa- 



tion 



Mo, = 



3eo/z 1 



F(T,,)Z(T,,) r(v)dv 



(4) 



where eo is the vacuum permittivity, jUgi is the permanent electric 
dipole moment, 5ui is the line strength, is the upper state 
energy, Z is the rotational partition function, and /rdv is the 
integrated optical thickness. We assumed a dipole moment of 
3.4 D for both N2H+ and N2D+ JHavenith et al. 19901 ). For the 
rotational transition — > /u - 1 of a linear molecule (like N2H^), 

For the N2H^ lines the optical thicknesses were derived from 
the Gaussian fits to the hyperfine components, and thus the 

integral J rdv can be replaced by ^^^ Avtq. Here Av is the 
linewidth of an individual hyperfine component, and tq is the 
sum of the peak optical thicknesses of all the seven components. 

The N2D^ column densities were calculated in two different 
ways: 1) as in the case of N2H^, and 2) using Eq. (1) with the 
approximation of optically thin line (r <K 1): 



(5) 



Again, we assumed that 77 - rjMB- The integrated opacity was 
estimated from the integrated Tmb of the main hyperfine group 
(54.3% of the total integrated intensity): 



rdv 



/ rMfidv 



(6) 



A comparison of the column density determination via the 
two methods shows that A^(N2D+), when using the first method, 
is 1.3 (IRAS 05405-0117), 0.7 (Ori B9 E), and 0.6 (Ori B9 N) 
times the value obtained using the second one. As the N2D^ line 
areas are somewhat uncertain, the N2D+ column densities deter- 
mined by using the first method have been adopted in this paper. 

The fractional N2H^ and N2D^ abundances, .«(N2H^) and 
ji:(N2D^), were calculated by dividing the corresponding column 
densities by A^(H2) from the dust continuum. For x(N2H+) the 
dust map was smoothed to 26'.'4, the resolution of the N2H^ 
observations. No smoothing was done in the case of N2D^, as 
the resolutions of the N2D^ and dust continuum observations 



are similar (16'.'0 and 18'.'6, respectively). The degree of deu- 
terium fractionation in N2H^ is defined as the column density 
ratio TJdeut = MN2D+)/A?(N2H+). 

The obtained H2 column densities are given in column (4) 
of Table |7] The N2H^ and N2D+ column densities, fractional 
abundances, and the values of /?deut are listed in Table [8] The 
uncertainties on A^(N2H^) and A^(N2D^) have been calculated 
by propagating the uncertainties on Tex, Ttot, and Av, and the 
uncertainties on A^(N2D^)/A^(N2H^) ratios are propagated from 
A^(N2H+)and A?(N2D+). 

4.4. Ionization degree and cosmic ray ionization rate 

The charge quasi-neutrality of plasma dictates that the number 
of positive and negative charges are equal. Since electrons are 
the dominant negative species, their fractional abundance nearly 
equals the sum of the abundances of positive ions, x(cations) ^ 
x(e). Thus, one may obtain a lower limit for the ionization frac- 
tion by summing the abundances of several molecular ions: 

x(e) > Jc(N2H+) + x(N2D+) + jc(H+) + jt:(H2D+) . (7) 

In the following we attemp to derive estimates for the cos- 
mic ray ionization rate and the fractional electron abundance us- 
ing the abundances of N2H^, N2D^, and H2D^ together with 
the reaction schemes and formulae presented in Crapsi et al. 
(2004) and CaselH et al. (2008). The rate coeflicients for the 
H3 + H2 isotopic system have been newly calculated by Hugo 
et al. (2009). We use these for the deuteration sequence H3 <-> 
H2D+ ^ D2H+ ^ D+ (see their Table VIII) . For other reac- 
tions the rate coefficients have been adopted from the UMIST 
database which is available at www . udfa . net. The main differ- 
ence between the Hugo et al. coefficients and those of Roberts et 
al. (2004) is that in the former, the effective backward rate coef- 

ficient, of the reaction Ht + HD ^ H2D^ + H2, and the cor- 

k-i 

responding coefficients for multiply deuterated forms of H3 are 
higher if the non-thermal ortho/para ratio of H2 (hereafter o/p- 
H2) is taken into account (Pagani et al. 1992 IGerhch et al. 20021 
Flower etal. 2 006a; Pag ani et al. 2009a[|Hugo et al. 2009| l. 

The ortho-H2D+ column density was derived towards Ori 
B9 E and N by Harju et al. (2006). Using their value, 
A^(o-H2D^) ~ 3.0 ■ 10'^ cm"^, and the H2 column densities 
derived here, we get x(o - H2D+) ^ 8.0 ■ 10""' and 4.1 ■ lO"'" 
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Table 8. N2H^ and N2D^ column densities, fractional abundances, and the column density ratio. 





A'(N2H+) 


A'(N2D+) 


x(NiH+) 


x(NiD+) 




Position 


[10'2 


[10" 


[10-'"] 


[10-"] 


Ri.ui = A'(N2D+)/MN2H+) 


IRAS 05405-0117 


9.14 ±0.08" 


3.19 ±0.37 


11.1 


4.9 


0.03 ± 0.004 


Ori B9 E 


4.54 ± 0.65" 


1.90 ± 1.39 


6.9 


5.0 


0.04 ± 0.03 


Ori B9 N* 


4.11 ± 1.51 


1.46 ±0.65 


3.9 


1.9 


0.04 ± 0.02 



" Harju et al. (2006) estimated slightly lower N2H+ column density, ~ 6 - 8 • 10'- cm--, toward IRAS 05405-0117 from the NzH+Cl - 0) 
data of Caselli & Myers (1994). 

' For the other velocity component A'(N2H+) = 1.56 ± 0.09 ■ 10'^ cm--, A'(N2D+) = 7.64 ± 5.44 • 10'" cm-^, and R^^^t = 0.05 ± 0.03. 



Table 9. Parameters derived in Sect. 4.4. 



Source 


-<H2D+) 


x(Ht) 




-^(e)„ 


<4e))l 




c, 




[10-'] 


[10-'] 


[10-'] 


[10-'] 


[10-'] 


[10-" s-'] 


[103 ^^-3/2 gl/2j 


Ori B9 E 


2.0 


1.8 


2.0/3.0 


6.4 


3.4 


2.0/1.0 


2.6 


Ori B9 N 


1.0 


0.9 


1.0/2.0 


6.3 


3.3 


2.5/1.3 


4.1 



" The first value is calculated from Eq. 0, whereas the second value is calculated from Eq. 
* This is the mean value between the lower and upper limit, where x(e); is calculated from Eq. 
The second value is derived by including HCO^ in the analysis (see text). 



towards Ori B9 E and N, respectively!! The ortho/para ra- 
tio of H2D^ (hereafter o/p-H2D^) depends heavily on 0/P-H2. 
According to the model of Walmsley et al. (2004, see their 
Fig. 3), the characteristic steady-state value of 0/P-H2 is ~ lO"'* 
in the density range n(H2) ~ 10^ - lO*" cm"^ appropriate for 
the objects of this study. This model deals with the situation of 
"complete depletion" and it is not clear how valid the quoted 
0/P-H2 is in less depleted gas. The recent results of Pagani et al. 
(2009a) suggest high values of 0/P-H2 (~ 4 ■ 10"^ - 6 ■ lO"^) in 
L183. 

For the moment we adopt the value 0/P-H2 = lO"'*. The 
effect of increasing this ratio will be examined briefly below. 
Assuming that o/p-H2D^ is mainly determined by nuclear spin 
changing collisions with ortho- and para-H2, the quoted 0/P-H2 
ratio implies an o/p-H2D^ of ~ 0.7 at T = 10 K. The total 
(ortho-Hpara) H2D^ abundances corresponding to this o/p ratio 
are x(H2D+) * 2.0 ■ lO"'' and 1.0 ■ 10"'^ towards Ori B9 E and N, 
respectively. 

The N2D^/N2H^ column density ratio which we denote by 
^deut, gives a rough estimate for the H2D^/H3 abundance ra- 
tio, denoted here by r. According to the relation /?deut ~ {r + 
2r^)/(3 + 2r + r^) derived by Crapsi et al. (2004; a more accurate 
formula is given in Eq. (13) of Caselli et al. (2008)), the values of 
/?deut given in Table[8]imply r ^ 0.08 for IRAS 05405-01 17, and 
r X! 0.1 1 for Ori B9 E and N. Using these r values we obtain the 
following fractional abundances ^(H^) - x(H2D^) / r » 1.8 ■ 10"^ 
(Ori B9 E) and a: 9.1 ■ 10"'^ (Ori B9 N). Substituting all the de- 
rived abundances into Eq. (|7]i we get the following lower limits 
for the degree of ionization: .^^(e) > 2.0 ■ lO"** in Ori B9 E, and 
> 1.0- 10"^ inOriB9N. 

Despite the fact that the clump associated with IRAS 05405- 
0117 does not stand out in the '^CO and C'^O maps, a mod- 
erate CO depletion factor of 3.6 near Ori B9 N has been de- 
rived ( [Caselli & Myers 1995tlCaseIIi et al. 2008) . This estimate 
is based on '■'CO(l - 0) observations made with the FCRAO 14- 
m telescope (HPBW 50"), and a total H2 column density, N{H2), 



' Caselh et al. (2008) derived Mo - H2D+) = 2.0/9.0 • 10'^ cm-^ 
toward position which is only 12'.' 7 southeast of our line observations 



derived from ammonia. By smoothing the LABOCA map to the 
resolution of 50", we obtain an average H2 column density of 
1.3 ■ 10^^ cm-^ around Ori B9 N. This is 2.6 times lower than 
the value adopted by Caselli et al. (2008). With this NiHi) the 
fractional CO abundance, x(CO), becomes 6.8 ■ lO"''. The corre- 
sponding CO depletion factor, /d, is only 1.4 with respect to the 
often adopted fractional abundance from Frerking et al. (19820 
The small value of Rdeut is consistent with a low CO depletion 
factor (e.g., Crapsi et al. 2004[ i 



In chemical equilibrium the fractional HI abundance is 



^H,/«(H2) + ^-lx(H2D+) 

Do 



(8) 



where is the cosmic ray ionization rate of H2, and 



Dq = ^ix(HD) + ^COJC(CO) + ^recOJc(e) + kgX(g ) + fcN2-«(N2) + ... 

The notation of Caselli et al. (2008) has been used here, i.e., ki 
and k-i are the forward and backward rate coefficients of the 
reaction mentioned above, and the other terms in Dq refer to the 
destruction of in reactions with neutral molecules (e.g., CO 
and N2) and in recombination with electrons and on negatively 
charged dust grains. 

By solving numerically Eqs. (8)-(10), and (13) of Caselli 
et al. (2008) together with our Eq. (O we obtain the follow- 
ing estimates for the fractional electron abundance and cosmic 
ray ionizations rate: x(e) - 6.4 ■ 10"', = 2.0 ■ lO""' s"' 
in Ori B9 E, and x(e) = 6.3 ■ lO"', (h. = 2.5 ■ lO"'^ s"' 
in Ori B9 N. Here we have used CO depletion factor 1.4, and 
the dust parameters {kg-, x{g^) quoted in Eqs. (11) and (12) of 
Caselli et al. which are based on a MRN dust grain size dis- 
tribution ( Ma this et al. 1977i i and effective grain recombination 
coefficients derived by Draine & Sutin (1987). The average num- 
ber densities, («(H2)), derived in Sect. 4.2. have been used for 
the cosmic ray ionization rates. The obtained values of ^h, are 



position Ori B9 N, assuming a critical density ric, 
respectively. 



10^ and 10^ cm- 



' The depletion factor is used in the text to express the fractional 
CO abundance with respect to the value 9.5 • 10-'. Adopting a higher 
reference abundance (see |Lacy et al. 1994| would not change the results 
of the calculations. 
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Velocity (km/s) 



Fig. 4. N2H+(1 - 0) spectra toward IRAS 05405-01 17 (top), Ori 
B9 E (middle), and Ori B9 N (bottom) after Harming smoothing. 
Hyperfine fits to the spectra are indicated with green lines. The 
residuals of the fits are shown below the spectra. Hyperfine fits 
to the other velocity component are indicated with red lines (see 
text). The small "absorption"-like feature at ~ 20 km s"' in the 
bottom panel is arfefact caused by frequency switching. 
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Fig. 5. N2D+(2 - 1) spectra toward IRAS 05405-01 17 (top), Ori 
B9 E (middle), and Ori B9 N (bottom) after Hanning smooth- 
ing. Hyperfine fits to the spectra are indicated with green lines. 
The lines under the spectra indicate the positions and relative 
intensities of the hyperfine components (see Table 2 in Gerin 
et al. (2001)). Undermost are plotted the residuals of the fits. 
Hyperfine fit to the other velocity component in the bottom panel 
is indicated with red line (see text). 



very similar to each other as is expected for such a nearby cores 
dWilliams et al. 19981 Bergin et al. 1999 and references therein). 
For comparison, Bergin et al. (1999) found that adopting 
- 5 ■ 10 s ' in their chemical model best reproduced their 
observations of massive cores in Orion. Note that the "standard" 
value often quoted in the literature is = 1-3 ■ 10"'^ s"'. Also 
the fractional electron abundances are clearly larger than those 
calculated from the standard relation x(e) ~ 1.3 ■ 10"^«(H2)"'^^ 
(cf. McKee 1989; McKee et al. 1993), where the electron frac- 
tion is due to cosmic ray ionization only and ^h, has its above 
mentioned standard value. The corresponding values would be 
Jt;(e) ~ 5 ■ 10 ^* (Ori B9 E) and ^e) ~ 3 ■ 10"^ (Ori B9 N). The 
mean value of the ionization degree found by Bergin et al. ( 1 999) 
for the massive cores in Orion is ~ 8 ■ 10"^. 

The parameters derived above depend on the adopted o/p-Ha 
which affects the backward rate coefficient A:_2, and A:_3 (see 



Caselli et al. 2008), and the CO depletion factor fo which affects 
the destruction of H3 . The fractional electron abundance can be 



decreased to ~ 8 ■ 
yields a ^h, of 1 .3 



10" by increasing 0/P-H2 to 2.4 ■ 10" (this 
10 s '). On the other hand, an increase in 
/d will lead to a higher x{e), but also to a lower ■ A solution 
where both x(e) and obtain the average values derived by 
Bergin et al. (1999) can be found by setting /d to 4.4 and 0/P-H2 
to 3.4 ■ 10^1 

However, the available observational data do not give 
grounds for abandoning the present bona fide /d value 1 .4. So 
the main uncertainty seems to be related to the unknown o/p- 
H2. Additional uncertainties to the ^h, values are caused by the 
rough density estimates, and by the fact that densities in the posi- 
tions observed in molecular lines are probably lower than the av- 
erage densities adopted in the analysis. The electron abundance 
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Fig. 6. SEDs of the sources SMM 3 (top), SMM 4 (middle), and 
IRAS 05045-0117 (bottom). 24 and 70 fim data points are de- 
rived from archival Spitzer/MIPS data, while the 870 fim mea- 
surement is performed with LABOCA. For IRAS 05045-0117 
we include also IRAS archival data (see Table |5]l. Icr error bars 
are indicated for Spitzer and LABOCA data points. The solid 
lines in all plots represent the sum of two (cold+warm) compo- 
nents (see columns (4) and (5) of Table|6]l. The "ripple" between 
- 30-60 yum in the SED of IRAS 05045-01 17 is due to simple 
logarithmic interpolation used to derive the luminosity. 



obtained assuming an o/p-H2 of 1.0 ■ 10 is likely to be an up- 



per limit. This o/p ratio corresponds to steady state in highly 
depleted dense gas with large abundances of and capable 
of efficient proton exchange with H2 (Flow er et al. 2 007). Their 
replacement by other ions in less extreme situations can sustain 
higher 0/P-H2 ratios. 

A substantial amount of CO implies the presence of HCO^ 
in the gas. By including the dissociative electron recombina- 
tion of HCO^, and the proton exchange reaction between NaH^ 
(or N2D^) and CO in the reaction scheme, the fractional HCO^ 
abundance can be solved. Through this estimate we get a slightly 
more stringent lower limit on the electron abundance than that 
imposed by Eq. (|7]i by demanding that 



x(e) > x(H+)(l + r) + x(N2H+)(l + Ra,ut) 

+ x(HCO+)(l+7?deut). 



(9) 



Here it has been assumed that N2H^ and HCO^ have similar de- 
grees of deuterium fractionation. By varying 0/P-H2 until elec- 
tron and the "known" cations are in balance we obtain with o/p- 
H2 = 2.7 10-^^ the lower Hmits x(e) > 3 ■ lO"** and x(e) > 2 ■ lO"** 
in Ori B9 E and N, respectively. The corresponding values of 



are 1.0- lO"'" s"' and 1.3 ■ 10" 



In both solutions HCO^ 



is more abundant than H^, whereas in the case of a large x{e) 
(small 0/P-H2) the HCO^ abundance lies between those of N2H^ 
and H2D+. 

The obtained values of the cosmic ray ionization rate vary 
smoothly with 0/P-H2, and all viable solutions point towards 
^H, ~ l-2-10-iS-i.InthemodelofMcKee(1989)these levels 
imply factional ionizations of ~ 1.1 - 1.6 ■ 10"^ at the density 10^ 
cm-^. These values lie between the lower and upper limits de- 
rived above. In what follows we assume that x{e) ~ 1 ■ 10"^, keep- 
ing in mind that true electron abundance is likely to be found 
within a factor of few from this value. 

There is also uncertainty about the most abundant ion. 
According to our calculation the electron abundance is an order 
of magnitude higher than the summed abundances of the positive 
ions H+, HCO+, and N2H+ for o/p - Hj = 1.0 ■ lO "*, whereas 
for the higher 0/P-H2 ratio ji:(HCO^) is comparable to .Ji:(e). This 
suggests that in the first case the reaction scheme misses the 
most abundant cation(s). In depleted regions with densities be- 
low 10^ cm^^ protons, H^, are likely to be the dominant ions 



( [Walmsley et al. 2004[ |Pagani et aklOOga). On the other hand, 
as discussed in Crapsi et al. (2004) and references therein, if 
atomic oxygen is abundant in the gas phase the major ion may 
be H30^. To our knowledge this ion has not yet been found in 
cold clouds (see also Caselli et al. 2008). When discussing the 
ambipolar diffusion timescale in Sect. 5.6 we will assume that 
the most abundant ion is either or HCO^ . 

We furthermore estimate the value of a constant, C,, that de- 
scribes the relative contributions of molecular ions and metal 
ions to the ionization balance (William s et al. 19981 their Eq. (4); 
Berginetal. 1999[ iPadoan et al. 2 004). The value of C, can be 



used to estimate the strength of the ion-neutral coupling in terms 
of the wave coupling parameter, W oc C, (see Sect. 5.7). In this 
analysis it is assumed that the electron abundance is determined 
by cosmic ray ionization balanced by recombination and it is ap- 
propriate for cores where Ay > 4 mag (i.e., ionization due to cos- 
mic rays dominates that resulting from UV radiation). Adopting 
the electron abundance 1 ■ 10"^ and ^h. 



ues of C, ^ 2.6 -4.1 ■ 10' cm 



-3/2 „1A2 



in our cores. These are 



similar to the value found by Bergin et al. (1999) for the mas- 
sive cores in Orion (3.6 ■ 10^ cm"-'^^ s'^^). McKee (1989) derives 
Ci = 3.2 ■ lO-' cm--'''^ s'^^ for an idealised model of cosmic ray 



12 



Miettinen et al.: N2H*, N2D* and dust emission in Ori B9 



ionization and Williams et al. (1998) obtained C, = 2.0 ■ 10^ 
(;jjj-3/2 gi/2 ^Qj. jow-mass cores. 

All the parameters derived in this Section are summarised in 
Tabled 



5. Discussion 

5.1. Nature ofsubmm sources in Ori B9 

By combining the submm LABOCA and far-infrared Spitzer 
data, we can distinguish starless cores from protostellar cores. 
In addition to the four IRAS sources in the region, two of the 
new submm sources, namely SMM 3 and SMM 4, are clearly 
associated with Spitzer point sources and are protostellar. The 
remaining six submm cores are starless. 

IRAS 05399-0121 was previously classified as a Class I pro- 
tostar (Bally et al. 2002 and references therein). However, tak- 
ing into account the rather low bolometric (18.5 K) and ki- 
netic temperatures (13.7 K, [Harju et al. 1993) , and high values 
of L,ubmm/iboi (2%) and M,ot/L-° (0.45 Mq/L^"), we suggest 
the source is in a transition phase from Class to Class I (see 
pontemps et al. 1996 [Froebrich 2005|. This source is associ- 
ated with the highly collimated jet HH 92 ( [Bally et al. 2002| l. 

The SED of IRAS 05413-0104 derived here is consis- 
tent with its previous classification as a Class object (e.g., 
[Cabrit et al. 2002 and references therein). The source is asso- 
ciated with the highly symmetric jet HH 212 dLee et al. 20()6l 
2007: ICodella et al. 20071 fSmith et al. 20071 ICabrit et al. 2007] l. 
IRAS 05412-0105 and IRAS 05405-0117 have very similar 
SEDs, and they, too, are likely to represent the Class stage. 
The weak line wings in the N2H^(1 - 0) hyperfine fines of IRAS 
05405-01 17 (see Fig.|4] top) could indicate the presence of out- 
flow from an embedded protostellar object. 

L,ubmm/Lbo\ ratio for both SMM 3 and SMM 4 is 11%. This 
together with low values of Tboi make these new submm sources 
Class candidates (e.g., ,Froebrich 2005 and references therein) 
that are deeply embedded in a massive, cold envelope. On a bolo- 
metric luminosity vs. temperature diagram these objects lie on 
the evolutionary track for a Class source with initially massive 
envelope (see Fig. 12 in Myers et al. (1998)). 

The stai-less cores SMM 1, 2, 5, 6, 7, and Ori B9 N, are likely 
to be prestellar as their densities are relatively high (0.6-5.8 ■ 10^ 



cm" ; see also Sect. 4.2). The 24 yum Spitzer source near SMM 
5 is probably not associated with this core. It lies rather far form 
the core centre and it is not detected at 70 fim (Fig.|2]l. 

There is an equal number of prestellar and protostellar cores 
in Ori B9. This situation is similar to that recently found by 
Enoch et al. (2008) in Perseus, Serpens, and Ophiuchus, and 
suggests that the lifetimes of prestellar and protostellar cores are 
comparable. Evolutionary timescales will be further discussed in 
Sect. 5.6. 



5.2. Mass distribution and core separations 

The spatial and mass distribution of cores are both important pa- 
rameters concerning the cloud fragmentation mechanism. Our 
core sample is, however, so small that it is not reasonable to 
study the properties of these distributions directly. Therefore, we 
only compared them with the distributions derived for another, 
larger core sample in Orion GMC by NW07. We make this com- 
parison particularly with Orion B North because the SCUBA 
850 yum map of Orion B North (see Fig. 2c in NW07) looks 
qualitatively similar to Ori B9. Orion B North also has deeper 
sensitivity and completeness limit than other regions studied by 



NW07, and besides it contains large number of cores. Fig. [T] 
presents the observed cumulative mass functions, which counts 
cores with mass less than M, i.e., N(M) - N(m < M)/Ntot, for 
both the core masses in Ori B9 and masses of prestellar cores in 
Orion B North derived by NW07. Note that the core mass func- 
tion (CMF) studied by NW07 is constructed by removing the 
Class I protostars from the sample, so that CMF includes only 
cores which have all their mass initially available for star forma- 
tion left. Correspondingly, we have excluded IRAS 05399-0121 
from our sample. We have also multiplied the core masses by 
required factors to compare to NW07 values, due to differences 
in assumed values of Tdust, k^, and distance (NW07 used the fol- 
lowing values: T^^^t - 20 K, k^^q = 0.1 m^ kg"', and d - 400 
pc). 
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Fig. 7. Normalised cumulative mass functions, N(M), for the 
prestellar cores in Ori B9 (solid line) and in Orion B North 
(dashed line) studied by Nutter & Ward-Thompson (2007). 

In order to determine if the two datasets are samples of the 
same core mass distribution, we carried out the Kolmogorov- 
Smirnov (K-S) test. The K-S test yields the maximum vertical 
difference between the cumulative distributions of D = 0.166, 
and probability of approximately 95% that the core mass distri- 
butions in Ori B9 and Orion B North are drawn from the same 
parent distribution. 

Fig. |8] (top) shows the observed core separation distribution 
and the distribution expected for the same number of randomly 
positioned cores over an identical area (0.22 deg-). The mean 
and median of the core separations in Ori B9 are log(r/AU) = 
5.47 + 0.09 (2.9 ± 0.6 ■ 10^ AU) and 5.42 (2.6 ■ 10^ AU), respec- 
tively. The quoted error for the mean correspond to the standard 
deviation. These values are similar to those of randomly posi- 
tioned cores, for which the mean and median are log(r/AU) = 
5.59+0.05 and 5.59±0.06, respectively. The quoted uncertainties 
are the standard deviation of the sampling functions. For core 
separation distribution in Orion B North studied by NW07, the 
corresponding values are 5.67+0.03 and 5.63, suggesting that the 
fragmentation scale is similar in both Ori B9 and Orion B North. 
Similar fragmentation scales and the fact that CMFs have re- 
semblance to the stellar IMF dGoodwin et al. 2008] l suggest that 
the origin of cores in these two regions is probably determined 
by turbulent fragmentation (e.g., Mac Low & Klessen 2004t 
Ballesteros-Paredes et al. 2007). The clustered mode of star for- 
mation in these two regions suggests that turbulence is driven 
on large scales (e.g., IKlessen 20011 ). Recently, Enoch et al. 
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(2007) found the median separations of log(r/AU) - 3.79, 
4.41, and 4.36 in nearby molecular clouds Ophiuchus, Perseus, 
and Serpens, respectively. The spatial resolution of the Bolocam 
(31") used by Enoch et al. at the distance of Ophiuchus, Perseus, 
and Serpens, is 0.02, 0.04, and 0.04 pc. The latter two are simi- 
lar to our resolution. The results suggest that the fragmentation 
scales in Perseus and Serpens are different from that in Orion. 

Fig. |8](bottom) shows the comparison between the observed 
nearest neighbour distribution and the distribution for randomly 
positioned cores. The mean and median of the nearest neighbour 
distribution in Ori B9 are log(r/ AU) = 4.75 ± 0.09 (5 .6 ± 1 .3 ■ 10* 
AU) and 4.62 (4.2 ■ lO'* AU), respectively. These values are 
rather different from those expected from random distributions, 
for which the mean and median are log(r/AU) = 5.08 ± 0.08 and 
5 .07 + 0. 1 1 , respectively. For core positions in Orion B North the 
mean and median are log(r/AU) = 4.46 + 0.03 and 4.35, respec- 
tively (NW07). Also this comparison supports the idea that the 
scale of fragmentation, and the amount of clustering are similar 
in Ori B9 and Orion B North. Note that the minimum observable 
separation is the beam size, i.e 18'.'6 or ~ 8.3 ■ lO"* AU at 450 pc. 
We also note that the source sample is too small to measure the 
significance of the clustering in Ori B9 based on the two-point 
correlation function. 



5.3. Sizes, shapes, and density structures of ttie cores 

Starless cores in Ori B9, for which the mean value of the 
deconvolved angular size in units of the beam FWHM is 
(Os/Obeam) - 2.5 + 0.8, are larger on average than protostellar 
cores ((^s/^eam) - 1 .6 + 0.2). Thcsc sizcs are similar to thosc re- 
cently found by Enoch et al. (2008) in Perseus ((fls/^beam) - 2.2 
and 1.6 for starless and protostellar cores, respectively). The 
mean axis ratios at half-maximum contours of starless and pro- 
tostellar cores are also 2.5 and 1.6, respectively (see Table |2] 
column (7)). This indicates that starless cores in Ori B9 are 
also more elongated on average than protostellar cores (cf. 
lOffner & Krumholz 20091 

The values of Os/Obeam can be used to infer the steep- 
ness of the core radial density profile (Youn g et al. 2003] 
[Enoch et al. 20081 l. According to the correlation between 
Os/Obeam and density -power-law index, p, found by Young et 
al. (2003, see their Fig. 27) a mean ^s/^beam values of 2.5 for 
starless and 1.6 for protostellar cores imply an average index of 
p ~ 0.9 - 1.0 and ~ 1.4-1.5, respectively. Moreover, Fig. 25 
of Young et al. (2003) suggest power-law indices < 1 for star- 
less cores and 1.1 - 1.6 for protostellar cores, consistent with 
those inferred by the average deconvolved angular source sizes. 
The low values of p ~ 1.0 for starless cores suggest that they 
are best modelled with shallower density profiles than the pro- 
tostellar cores. These results are in agreement with those found 
by Ward-Thompson et al. (1999) using 1.3 mm dust continuum 
data, and CaselH et al. (2002b) using N2H+(1 - 0) maps. 



5.4. Deuterium fractionation and depletion in the IRAS 
05405-0117 region 

The N2D^/N2H^ column density ratio, Rdeuu is supposed 
to increase strongly as the core evolves (Caselli 2002( 
Crapsi et al. 2005a| thei r Fig. 5; iFontani et al. 2006; 
Emprechtinger et al. 2009} but see [Roberts & Millar 20070 . 
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This can be understood so that the abundances of H^, and its 
deuterated forms which transfer deuterium to other molecules, 
increase with increasing density due to molecular depletion and 



log(r/AU) 

Fig. 8. Top: Observed core separation distribution (solid line) 
compared with the expected distribution for random distribution 
of the same number of sources as the observed sample over an 
identical area (dashed line). Bottom: Observed nearest neigh- 
bour distribution (solid line) compared with the expected distri- 
bution for random distribution of the same number of sources as 
the observed sample over the same area (dashed hne). 



a lower degree of ionization. Crapsi et al. (2005a) suggested 
that prestellar cores are characterised with Rdeut > 0.1, whereas 
starless cores with Rdem < 0.1 are not necessarily so dense 
that CO would be heavily depleted dRoberts & Millar 2007"] l. It 
should be noted that in cores without internal heating sources 
the degree of deuterium fractionation is likely to increase 
inwards as the density increases and temperature decreases due 
to attenuation of starlight. This temporal and radial tendency 
is likely to be reversed during the core collapse because of 
compressional heating and the formation of a protostar (e.g., 
lAikawa et al. 2008 a. see also Fig. 3 in Emprechtinger et al. 
(2009)). 

The positions studied here have Rdeut ~ 0.03-0.04. This is ~ 
2-3 10^ times larger than the cosmic D/H elemental abundance 
of ~ 1.5 • 10-5 ( [Linsky et al. T995j 2006: lOliveira et al. 2003b . 
Also, the H2D^/H3 ratios we derived are ~ 7 ■ 10^ times larger 
than the cosmic D/H ratio. Our Rdem values are similar to those 
found by Crapsi et al. (2005a) toward several low-mass starless 
cores, and to those found by Emprechtinger et al. (2009) to- 
ward Class sources. Like Emprechtinger et al. (2009), we find 
that the deuterium fractionation of N2H^ in protostellar cores. 
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which takes place in the cold extended envelope, is similar to 
that in prestellar cores. It has been found that the values of Rdeut 
toward high-mass star-forming cores are (usually) lower than 
those found in the present study ( IFontani et al. 20061 see also 
Emprechtinger et al. 2009[ l. This conforms with the fact that our 



sources are low- to intermediate mass star-forming cores. 

As discussed by Walmsley et al. (2004) and Flower et al. 
(2006a), the H2D^ abundance depends (inversely) on the or- 



H+-HHD 



tho:para ratio of H2, because the reaction HaD^ -1- H2 
is rapid between ortho forms. The ortho:para ratio of H2 de 
creases with time and gas density, and is therefore large at early 
stages of core evolution. Consequently, a relatively high degree 
of deuterium fractionation is a sign of matured chemistry char- 
acterised by a low ortho:para ratio of H2 and probably a high 
degree of molecular depletion. Low CO depletion factor of 1 .4 
close to N2H^ peak Ori B9 N (see Sect. 4.4) is consistent with 
the 7?deut value of Ori B9 N (e.g., |Crapsi et al. 2004[ see also 
Fig. 4 in Emprechtinger et al. (2009)). The ortho-H2D^ detec- 
tion towards Ori B9 E and N suggests an evolved chemical stage 
and tells of a longlasting prestellar phase. The non-detection to- 
ward IRAS 05405-01 17 can be explained by a lower ortho-H2D+ 
abundance due to the central heating by the protostar 

5.5. Evidence for a N2H* "hole" and chemical differentiation 

The N2H+ map of Caselli & Myers (1994, see their Fig. 2) and 
submm dust continuum map of the clump associated with IRAS 
05405-01 17 (Fig. [3]) are not very much alike. The strongest dust 
continuum peak, SMM 4, does not stand out in N2H^. Moreover, 
the northern N2H^ maximum, Ori B9 N, seem to be shifted with 
respect to the northern dust peak, and the N2H^ peak Ori B9 E 
does not correspond to any dust emission peak. 

To determine whether the ~ 1 .5 times higher resolution of 
the LABOCA 870 fim map relative to the N2H+ map of Caselli 
& Myers (1994) contributed to the different appearance of the 
dust continuum and N2H^ maxima, we smoothed the LABOCA 
map to a resolution similar to that of the N2H^ map (27"). The 
smoothed 870 fim map, however, still shows the same differ- 
ences between the dust continuum and N2H^. 

The Class candidate SMM 4 (see Sect. 4.1 and 5.1) can 
represent an extreme case of depletion where also N2H^ has 
disappeared from the gas phase due to freeze out on to the 
dust grain surfaces. There is some previous evidence for N2H^ 
depletion in the centres of chemically evolved cores, such as 
B68 (|Bergin et al. 2002| l, L1544 tCaselli et al. 2002ai L1512 



dLee et al. 2003D , and Li521F ( [Crapsi etal."20 04). Example of 
the N2H+ depletion toward Class source is IRAM 04 1 9 1 + 1 522 
in Taurus (Be lloche & Andre 200 4). Pagani et al. (2005) found 
clear signs of moderate N2H^ depletion in the prestellar core 
L183 (see also [Pagani et al. 2007|). Also, Schnee et al. (2007) 



found clear evidence of N2H'^ depletion toward the dust centre 
of TMC-IC. Note that SMM 4 is probably not warm enough 
for CO to evaporate from the grain mantles (~ 20 K, e.g., 
lAikawa et al. 2008al l. so it is unlikely that CO, which is the main 
destroyer of N2H+ (through reaction N2H++CO ^ HCO++N2), 
would have led to disappearance of N2H^ from the gas phase. 

To study the chemical differentiation within the clump, we 
compare our previously determined NH3 column densities with 
the present N2H^ column densities. The integrated NH3 (1, 1) 
and (2, 2) intensity maps of the clump (see Appendix A in Harju 
et al. (1993)) show roughly the same morphology as the submm 
map. The NH3 column densities toward IRAS 05405-0117, Ori 
B9 E, and Ori B9 N are 11.9 ± 1.3 ■ lO'"*, 7.2 + 1.3 ■ 10'^, and 



9.7 ±4.6- lO'** cm"^, respectively. The corresponding NH3/N2H^ 
column density ratios are about 130+ 14, 159 ±37, and 236 ±87. 
These values suggest that NH3/N2H+ abundance ratio is higher 
towards starless condensations than towards the IRAS source. 
Hotzel et al. (2004) found a similar tendency in B217 and L1262: 
the NH3/N2H^ abundance ratios are at least twice as large in the 
dense starless parts of the cores than in the regions closer to the 
YSO (see Casell i et al. 2002bl for other low-mass star-forming 
regions). The same trend is also found in the high-mass star- 
forming region IRAS 20293H-3952 (IPalau et al. 20071). This is in 
accordance with chemistry models (lAikawa et al. 20051 ) and pre- 
vious observations (ITafalla et al. 20041 which suggest that NH3 
develops slightly later than N2H^, and can resist depletion up to 
higher densities. It should be noted that models by Aikawa et al. 
(2005) reproduce the observed enhancement of the NH3/N2H^ 
ratio by adopting the branching ratio for the dissociative re- 
combination of N2H^ as measured by Geppert et al. (2004; i.e., 
N2H^ + e -> NH + N accounts for 64% of the total reaction). 
However, this branching ratio has been retreated by the same au- 
thorS and thus it is not clear at the moment what is actually 
causing the increase of NH3/N2H+ abundance ratio. 

5.6. Core evolution: quasi-static vs. dynamic 

The degree of ionization in dense cores determines the impor- 
tance of magnetic fields in the core dynamics. The ionization 
fractions in low-mass cores are found to be 10"^ < x(e) < 10"^ 
dCaselfi et al. 1998IIWiniams et al. 199 8'). The physical origin of 
the large variations in x(e) is not well understood, though varia- 
tions in or appropriate values of metal depletion are assumed 
jPadoan et al. 2004D . Padoan et al. (2004) suggested that the ob- 
served variations in x(e) can be understood as the combined ef- 
fect of variations in core age, extinction, and density. 

Fractional ionizations can be transformed into estimates of 
the ambipolar diffusion (AD) timescale, tad- We have used 
Eq. (5) of Walmsley et al. (2004) for this purpose. Assuming 
that is the dominant ion (see Sect. 4.4), one obtains t/^j) 8 ■ 
10'^x(e)yr. If the dominant ion is HCO+, tad ~ 1.3- 10''*x(e) yr, 
i.e., ~ 60% longer than in the former case. Using the elec- 
tron abundance x(e) = 1 ■ 10"^ we obtain that tad ~ 10^ 
yr. This timescale is roughly 70 and 100 times longer than 

the free-fall time (rff ~ 3.7 ■ 10^ (n(H2)[cm-3])"'^^ yr) of Ori 
B9 E and N, respectively. Since tad > t^, the cores may 
be supported against gravitational collapse by magnetic fields 
and ion-neutral coupling. The magnetic field that is needed 
to support the cores can be estimated using the relation be- 
tween the critical mass required for collapse and the magnetic 
flux (see Eq. (2) in Mouschovias & Spitzer (1976)). Using the 
masses and radii from Table |2l we obtain a critical magnetic 
field strength of ~ 80 juG for Ori B9 E/IRAS 05405-0117 and 
~ 200 juG for Ori B9 N. These are rather high values com- 
pared to those that have been observed (Troland et al. 1996t 
Crutcher 1999 , Crutcher & Troland 2000, Crutc her et al. 200^ 
iTurner & Heiles 20061 ITroland & Crutcher 20081) . 

According to the "standard" model of low-mass star 
formation, tap / Tff ~ 10 (see, e.g., IShu et al. 19871 
ICiolek & Basu 20011 and references therein). Since AD is 
generally a slow process, the core evolution toward star for- 
mation occur quasi-staticly. The chemical abundances found 
in the present study (x(N2H+) ~ lO""', x(NH3) ~ 10"^ 



Their recent laboratory experiment suggest that the above men- 
tioned branching ratio is only 10% (see lAikawa et al. 2008bt . 
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see Sect. 5.8) are consistent with chemical 
dynamically young, but chemically evolved 
(|Bergin & Langer 1997 



yr) source 

lAikawa et al. 2005 ; Shirley et al. 2005fsee also Ki rk et al. 20071 



models for a 
(age > 10^ 
[Roberts et al. 20041 



and references therein). This supports the idea that the sources 
have been static or slowly contracting for more than 10^ yr, and 
conforms with the estimated AD timescales. 

On the other hand, the equal numbers of prestellar and 
protostellar cores suggest that the prestellar core lifetime 
should be similar to the lifetime of embedded protostars. 
Since the duration of the protostellar stage is ~ few • 
10^ yr (e.g., [Ward-Thompson et al. 2007} IHaj chelL et al. 20071 
IGalvan-Madrid et al. 20071 Enoch et al. 20081) . the prestellar 
core evolution should be rather dynamic and last for only a 
few free-fall times, as is the case in star formation driven 
by supersonic turbulence (e .g., IMac Low & Klessen 20041 
IBallesteros-Paredes et al. 2007l l. This seems to contradict with 
the above results of AD timescales. However, in order to recog- 
nise the cores in the submm map, they are presumed to be in the 
high-density stage of their evolution. Thus, the short statistical 
lifetime deduced above is still consistent with the quasi-static 
evolution driven by AD, if we are only observing the densest 
stages of a longer scale core evolution (e.g., Enoch et al. 20081 
ICrutcher et al. 2009t . Also, the dynamic phase in the core evo- 
lution with Tad being only a few free-fall times might be appro- 
priate for magnetically near-critical (the mass-to-magnetic flux 
ratio being ~ 80% of the critical value) or already slightly super- 
critical cores when rapid collapse ensues jCiolek & Basu 20011 
see also lTassis & Mouschovias 20041 ). 

5. 7. Linewidths and turbulence 

The N2D^ linewidths in Ori B9 E and N are significantly 
narrower than the N2H^ linewidths (by factors of ~ 1.5 and 
~ 1.9, respectively). Crapsi et al. (2005a) found similar trend 
in several low-mass starless cores (see their Table 4). This 
is probably due to the fact that NaD^ traces the high den- 
sity nuclei of starless cores, where non-thermal turbulent mo- 
tions are expected to be insignificant (e.g., Andre et al. 2007, 
Ward-Thompson et al. 2007| l. 

The non-thermal component dominate the N2H^ linewidths 
in the observed positions (the thermal linewidth of N2H^ is about 
0.126 km s"' at 10 K, and thus Avnt/Avt ~ 2). However, the 
level of internal turbulence, as estimated from the ratio between 
the non-thermal velocity dispersion and the isothermal speed of 
sound (e.g. .' Kirk et al. 20071 . is not dynamically significant. 

Using Eq. (7) of Williams et al (1998) and the derived values 
for the molecular/metal ion-contribution constant C, and cosmic 
ray ionization rate (see Table |9]l, we see that non-thermal 
NaH^ line broadening in the observed positions can be explained 
in part by magnetohydrodynamic (MHD) wave propagation. The 
large wave coupling parameter in our sources {W » 1), sug- 
gest that the coupling between the field and gas is strong and the 
waves are not suppressed. The derived values of W (~ 30-60, in 
the case of minimum turbulence) are in agreement with TAD/^ft 
ratios (I Williams et al. 1998] l. Also, the estimated degrees of cou- 
pling between the magnetic field and gas conforms with the sus- 
ceptibility to fragmentation (Bergin et al. 1999). 



Caselli & Myers (1995) analysed ammonia cores in the 
Orion B GMC and found an inverse relationship between core 
linewidth and distance to the nearest stellar cluster The nearest 
stellar cluster to Ori B9 is NGC 2024 at the projected distance of 
5.2 pc (see Sect. 1.1.), so its role as driving external turbulence 
to the region is probably not significant. 



5.8. Formation of a small stellar group In Ori B9 

Internal turbulence or gravitational motions in the massive 
molecular cloud core may promote fragmentation of the 
medium. This can easily generate sheets and filaments (e.g.. 
Casein & Myers 19 95; And re et al. 20081 ). The collapse of these 
elongated clumps most probably results in the formation of a 
small stellar group or a binary system rather than a single star 
(e.g., ILaunhardt et al. 19961 ). Only the densest parts of the fila- 
ments, the dense cores, are directly involved in star formation. 
It is unclear at the present time whether the collapse of an indi- 
vidual prestellar core typically produces single stars or multiple 
protostellar systems (see I Andre et al. 2 008). 

The total mass of gas and dust of the clump associated with 
IRAS 05405-01 17 as derived from the dust continuum emission 
is ~ 14 Mq, and it has elongated structure with multiple cores 
(local maxima in the filament are separated by more than one 
beam size, see Fig. [3]). The previous mass estimates by Harju et 
al. (1993) based on NH3 were much higher: ~ 50 Mq derived 
from A^(NH3) distribution, and ~ 310 Mq derived from peak 
local density. The uncertainty in the abundance[]3 and lower res- 
olution used are certainly affecting the estimation of the mass 
from NH3. However, the clump has enough mass to form a small 
stellar group. 

The kinetic temperature, velocity dispersion and the frac- 
tional H2D^ abundance in the clump are similar to those in 
the well-studied prestellar cores, e.g., L1544 and L183, where 
strong emission of H2D+ line has been detected previously (see 
Harju et al. 2006| and references therein). The masses, sizes, rel- 
atively high degree of deuteration and the line parameters of 
the condensations indicate that they are low- to intermediate- 
mass dense cores (cf. Fontani et al. 20081 ). IRAS 05405-0117 
and SMM 4 are likely to represent Class protostellar cores (see 
Sect. 5.1), whereas the subsidiary cores, e.g., Ori B9 N, are in an 
earlier, prestellar phase. 

6. Summary and conclusions 

We mapped the Ori B9 cloud in the 870 /zm dust continuum 
emission with the APEX telescope. We also observed N2H^(1 - 
0) and N2D+(2 - 1) spectral line emission towards selected po- 
sitions in Ori B9 with the IRAM 30 m telecope. These obser- 
vations were used together with archival Spitzer/MIPS data to 
derive the physical characteristics of the cores in Ori B9 and the 
degree of deuterium fractionation and ionization degree within 
the IRAS 05405-0117 clump region. The main results of this 
work are: 

1. The LABOCA field contains 12 compact submm sources. 
Four of them are previously known IRAS sources, and eight of 
them are new submm sources. All the IRAS sources and two of 
the new submm sources are associated with the Spitzer 24 and 
70 fim sources. The previously unknown sources, SMM 3 and 
SMM 4, are promising Class candidates based on their SEDs 
between 24 and 870 yum. There is equal number of starless and 
protostellar cores in the cloud. We suggest that the majority of 
our starless cores are likely to be prestellar because of their high 
densities. 

2. The total mass of the cloud as estimated from the 2MASS 
near-infrared extinction map is 1400 Mq. The submm cores con- 
stitute about 3.6% of the total cloud mass. This percentage is in 



" Harju et al. assumed that x(NH3) ~ 3 • 10"**. Using the H2 column 
densities from the dust continuum we derive the values of x(NH3) ~ 
1 - 2 • 10"^ in our line observation positions. 
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agreement with the observed low values of star formation effi- 
ciency in nearby molecular clouds. 

3. Mass distribution of the cores in Ori B9 and in Orion B 
North studied by Nutter & Ward-Thompson (2007) very likely 
represent the subsamples of the same parent distribution. The 
CMF for the Orion B North is well-matched to the stellar IMF 
dGoodwin et al. 20^ 8^. Also the core separations in these two re- 
gions are similar, indicating that the fragmentation length scale 
is similar Since the fragmentation length scales are alike, and 
the CMFs have resemblance to the IMF, the origin of cores could 
be explained in terms of turbulent fragmentation. The clustered 
mode of star formation in these two different regions suggest that 
turbulence is driven on large scales. 

4. On average, the starless cores are larger and more elon- 
gated than the protostellar cores in Ori B9. The observed mean 
angular sizes and axis ratios suggest average density-power-law 
indices p ~ I and ~ 1 .5 for starless and protostellar cores, re- 
spectively. 

5. The fractional NaH^ and N2D^ abundances within the 
clump associated with IRAS 05405-0117 are ~ 4 - 11 ■ 10""' 
and 2 - 5 ■ lO"", respectively. The A?(N2D+)/A^(N2H+) column 
density ratio varies between 0.03-0.04. This is a typical degree 
of deuteration in low-mass dense cores and conform with the 
earlier detection of HaD^. There is evidence for a N2H^ "hole" 
in the protostellar Class candidate SMM 4. The envelope of 
SMM 4 probably represents an extreme case of depletion where 
also N2H^ has disappeared from the gas phase. 

6. The ionization fraction (electron abundance) in the posi- 
tions studied is estimated to be x(e) ~ 10"^. There is uncertainty 
about the most abundant ionic species. The most likely candi- 
dates are and HCO^. The cosmic ray ionization rate in the 
observed positions was found to be ~ 1 - 2 ■ 10"'* s"' . 

7. There seems to be a discrepancy between the chemical age 
derived near IRAS 05405-0117 and the statistical age deduced 
from the numbers of starless and protostellar cores which sug- 
gest that the duration of the prestellar phase of core evolution is 
comparable to the free-fall time. The statistical age estimate is, 
however, likely to be biased by the fact that the cores detected in 
this survey are rather dense (n(H2) > 10^ cm"^) and thus repre- 
sent the most advanced stages. 
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